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Abstract: We calculate resummed perturbative predictions for gaps-between-jets 
processes and compare to HERA data. Our calculation of this non-global observ- 
able needs to include the effects of primary gluon emission (global logarithms) and 
secondary gluon emission (non-global logarithms) to be correct at the leading loga- 
rithm (LL) level. We include primary emission by calculating anomalous dimension 
matrices for the geometry of the specific event definitions and estimate the effect of 
non-global logarithms in the large A*";, limit. The resulting predictions for energy flow 
observables are consistent with experimental data. 
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1. Introduction 

The subject of interjet energy flow [1] has attracted considerable interest ever since 
it was proposed [2, 3] as a way to study rapidity gap processes using the tools of 
perturbation theory. Rapidity gap processes are defined as processes containing two 
high pt jets with the region of rapidity between the jets containing nothing more than 
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soft radiation. This region is known generically as the gap. The presence of a range 
of scales offers a chance to study the interface between the soft, non-perturbative 
scales and the hard, perturbative scales of QCD. 

In this paper we will calculate the perturbative contribution to gaps-between- 
jets cross sections in photoproduction at HERA, which have been measured by the 
ZEUS [4, 5] and HI [6] collaborations. A feature of the recent analyses is the use 
of a clustering algorithm to define the hadronic final state and hence the gap. The 
restriction of transverse radiation in a region of phase space, defined as fl and directed 
away from the observed jets and the beam directions, produces logarithms at each 
order of QCD perturbation theory of the intcrjct energy flow, Qq, over some hard 
scale, Q. The precise deflnition of the restricted region, or gap, is totally free and 
in this work we are interested in the gap region defined by experimental rapidity 
gap analyses. The source of the large logarithms is twofold. The so-called primary 
(or global) logarithms arise from radiation emitted directly into Q; these wide-angle 
gluons decouple from the dynamics of the colour-singlet jets and are described by an 
effective, eikonal theory [7-11]. The second source of leading logarithms arise from 
gluons emitted outside of the gap region, an area of phase space generically denoted 
as n, which subsequently radiate into Q. These terms are known as non-global 
(secondary) logarithms, or NGLs [12-15]. 

The primary logarithms arc rcsummcd using the formalism of Collins, Soper and 
Sterman (CSS) [7, 16-18]. In this method the cross section is factorised into a soft 
part describing the emission of soft, wide angle gluons up to scale Qq and a hard 
part, describing harder quanta. A unique feature of QCD is that the soft and the 
hard functions are expressed as matrices in the space of possible colour flow of the 
system. The scale invariance and factorisation properties of the cross section are 
then exploited to resum primary logarithms of Qn/Q- This resummation is driven 
by the ultraviolet pole parts of eikonal Feynman graphs and we write the resummed 
cross section in terms of the eigenvalues of fi-dependent soft anomalous dimension 
matrices. These matrices are known for gap deflnitions based on the cone deflnition 
of the flnal state [9, 10] and for a gap deflned as a square patch in rapidity and 
azimuthal angle [11]; here we are interested in gaps deflned in terms of the clustering 
algorithms employed in the recent analyses. Hence we are required to calculate the 
corresponding anomalous dimension matrices. 

The NGLs [12, 13] are unable to be incorporated into the resummation of the 
primary logarithms, because the gluon emission patterns that produce the NGLs 
are sensitive to underlying colour flows not included in the formalism. The effect 
of NGLs, which is a suppressive effect, on energy ffow processes has been studied 
using numerical methods in the large N,. limit and overall factors describing their 
effect have been extracted for a two jet system, both without [13] and with [14] 
clustering. This factor is not directly applicable to the 4 jet systems'*" relevant in the 
^Note that for a two- to- two process the incoming and outgoing partons radiate, so we consider 
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photoproduction of jets but, in the lack of a four jet formalism, we nevertheless use 
the two-jet factor in our predictions. 

Our aim is to derive LL resummed predictions for the gap cross section, with 
primary logarithms correct to all orders and secondary logarithms correct in the 
large Nc limit. The gap cross section will follow the HERA analyses and demand two 
hard jets, defined using the kt clustering algorithm [19-21], and we will closely follow 
the HI and ZEUS gap definition. The technical aspects of soft gluon resummation 
give a strong dependence on the gluon emission phase space, and hence a considerable 
part of our work will be concerned with the calculation of soft gluon effects for the 
specific detector geometry of the HI and ZEUS experiments. 

The organisation of this paper is as follows. Section ^ describes, in detail, the 
energy flow analyses of HI and ZEUS. We describe the experimental cuts employed 
and the range of measured observables. We also discuss the theoretical implementa- 
tion of the inclusive kt algorithm employed to define the hadronic final state and the 
impact on soft gluon resummation. Section |^ describes the theoretical definition of 
our cross section and we employ the standard QCD factorisation theorems to write 
it as the convolution of non-perturbative parton distributions and a short-distance 
hard scattering function. We then proceed to refactorise the hard scattering function 
and exploit this factorisation to resum the large interjet logarithms. Section ^ then 
derives the soft anomalous dimension matrices for the kt defined final state and in 
section || we present detailed predictions of rapidity gap processes and compare to 
the HI data. Finally we draw our conclusion in section ^. We find that our descrip- 
tion of the data is good, although the approximate treatment of NGLs results in a 
relatively large normalisation uncertainty. 



2. The HERA energy flow analyses 

In this section we will outline the experimental analyses of the photoproduction 
of gaps-between-jets processes and discuss the experimental cuts and rapidity gap 
observables. We will also describe the clustering algorithm used to define the final 
hadronic state in the more recent ZEUS [5] and HI [6] analyses. 

The data for these events were collected when HERA collided 27.6 GeV positrons^ 
with 820 GeV protons, giving a centre of mass energy of ^/s ~ 300 GeV. Following 
the jet-finding phase, which we will comment on later, the total transverse energy 
flow between the two highest Et jets, denoted E^^^ , is calculated by summing the 
transverse energy of all particles that are not part of the dijets in the pseudorapidity 
region between the two highest jets. An event is then defined as a gap event if the 
energy is less than some energy cut E^^"^ = Qq. A gap fraction is then calculated 

the process to be of "four jet" type, although only two jets are seen in the final state. 
^The positron energy varied a negligible amount between the two sets of analyses. 
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by dividing the cross section at fixed E^^ by the inclusive cross section. The ZEUS 
collaboration performed a rapidity gap analysis several years ago [4] using the cone 
algorithm for the jet definition and presented the gap fraction at Qq = 0.3 GeV. We 
consider this value of Qn as being too small for our perturbative analysis and will 
not make any predictions for this data set. The more recent HI and ZEUS analyses 
used the kt definition of the final state and both collaborations presented the gap 
fraction at four different values of Q^, as shown in table We will make predictions 
and compare to data for the HI data sets and, due to the fact that the ZEUS data is 
still preliminary, confine ourselves to making predictions for the ZEUS analysis. We 
have summarised the cuts used in table [^. 





HI 


ZEUS 




> 6.0 GeV 


> 6.0 GeV 


Ef' 


> 5.0 GeV 


> 5.0 GeV 


^jetl 


< 2.65 


< 2.4 


^jet2 


< 2.65 


< 2.4 


Ar/ 


2.5 < Ar] <4: 


2 < Ar/ < 4 


mi 


N/A 


< 0.75 


y 


0.3 <y < 0.6 


0.2 < y < 0.85 




< 0.01 GeV2 


< 1 GeV2 


jet def. 


kt 


kt 


gap def. 


Ay = Ai] 


Ay = Af] 


R 


1.0 


1.0 


Qn 


0.5,1.0,1.5,2.0 GeV 


0.5,1.0,1.5,2.0 GeV 



Table 1: The experimental cuts used for the HERA analyses. 



2.1 The kt algorithm 

Of special interest to those going about soft gluon calculations is the method used 
to define the hadronic final state, the reason being that this jet-finding process de- 
termines the phase space for soft gluon emission; the method used in the HI and 
ZEUS data sets is the inclusive kt algorithm [19-21]. In this algorithm the final state 
is represented by a set of "protojets" i with momenta and works in an iterative 
way, grouping pairs of protojets together to form new ones. The aim is to group 
almost-parallel protojets together so that they are part of the same protojet. Once 
certain criteria are met, a protojet is considered a jet and is not considered further. 
Here we follow the so-called inclusive scheme used at HI and ZEUS which depends 
on the parameter R, normally set to unity. If we assume that any radiation into the 
gap is much softer than any parent radiation, then this radiation with E^ < -E^* will 
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be merged into the jet (with kinematical variables {fjjet, 0jet)) if it satisfies 

(r?r-r7jet)'+(0r-0jet)'<i?', (2.1) 

where we denote the kinematical variables of the radiated gluon by {r]r,(f>r)- Once 
merged, a gluon will be pulled out of the gap and can no longer produce a primary 
or secondary logarithm. The gap is defined as the interjet region minus the region 
of clustered radiation around the jets and may contain soft protojets. The gap 
transverse energy is then defined by the (scalar) sum of the protojets within the gap 
region, 771 < 77 < 772- 

The kt gap definition can be contrasted to the older ZEUS analysis [4] , which used 
the well known cone definition of the final state with R — 1.0. The gap transverse 
energy is then defined as the scalar sum of the hadrons within it, 771 + i? < 77 < 772 — -R- 

3. Factorisation, refactorisation and resummation of the cross 
section 

In this section we will exploit the standard factorisation theorems of QCD to write 
down the dijet production cross section from the interaction of a proton and a 
positron. We will then refactorisc the hard scattering function into the product 
of two matrices in the space of possible hard scattering colour fiow, one matrix de- 
scribing soft gluons radiated into the gap region and the other a hard scattering 
matrix. The renormalisation properties of the cross section are then used to resum 
primary interjet logarithms, and write the result in terms of the eigenvalues of the 
matrix of counterterms used to renormalise the soft function. In the following section 
we will calculate these matrices and their eigenvalues. 

3.1 Photoproduction cross sections 

The scattering of positrons and protons at HERA proceeds predominantly through 
the exchange of photons with very small virtuality and produces a large subset of 
events with jets of high transverse momentum, Et- The presence of this large scale 
allows the application of the perturbative methods of QCD to predict the cross section 
for multiple jet production. This process is otherwise known as jet photoproduction. 

The leading order (LO) QCD contribution can be divided into two types [9]. 
The first is the direct process in which the photon interacts directly with a parton 
from the proton and proceeds through either the Compton process, 7g gq, or 
the photon-gluon fusion process, 'yg qq. The second contribution is the resolved 
contribution, in which the virtual photon fluctuates into a hadronic state that acts 
as a source of partons, which then scatter off the partonic content of the proton. 
Therefore the reaction proceeds through standard QCD 2 — > 2 parton scattering 
processes. Note that the precise determination of the partonic content of the photon 
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is a very open question and there is a relatively large error associated with the 
photonic parton densities. The spectrum of virtual photons is approximated by the 
Weizacker-Williams [22] formula, 

« (i + (i-yn,„„/eax(i-?/) 



F,/e{y) = 1^^-^ ^ log ^^^^^^^^V^ ' (3-1) 



where nie is the electron mass, y is the fraction of the positron's energy that is 
transfered to the photon, and Q^ax maximum virtuality of the photon, which 

is determined by the experimental cuts employed in the analyses. Then, by using 
the equivalent photon approximation, the cross section for the process e'^p — > e'^X 
is given by the convolution 

dyF^/,{y)da{^p^X), (3.2) 

where we write d(T{'~fp X) for the cross section of 7p X. The centre of mass 
energy squared for the photon-proton system is W"^ = ys, where s is the centre of 
mass energy squared for the positron-proton system. At HERA, s ~ 90, 000 GeV^ 
and the values for y^^in and ymax are determined by the experimental analyses. We 
can now write down the specific expression for the production of two high Et jets 
from the photon-proton system, which is written as a sum of the direct and resolved 
contributions, 

da^+p{s, i, Arj, as{Hr), Qn) = / dy -F^/e(y) {da^p{s^p, i, Arj, asi^r), Qn) 

^ Vrain 

+da;;^(s^p, t, Ar/, as^iir). Qn)) , (3.3) 

where we denote the 4-momentum transfer squared in the hard scattering as t. We 
define the rapidity separation and difference of the two hard jets by 

A?7 = \■r]x-^2\, 

r]jJ-\{r]i^r]i). (3.4) 

At this point we can appeal to the collinear factorisation theorems of QCD and, by 
working in the 7p frame, write down factorised forms for the direct and resolved 
cross sections. The factorised direct cross section is 



d(7, 



dir 



JpJl,J2 



x—^{s,t,Ar),as{iir),Qn,l^f), (3.5) 
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and the factorised resolved cross section is 



x-^{s,t,Ar),as{Hr),Qn,l^f), (3.6) 

which are written in terms of the jet rapidity, rj, in the partonic centre-of-mass 
frame, and we write the factorisation scale and the renormalisation scale as /// and 
fjLr respectively. Note that fj — A77/2, s — XpW^ for the direct case and s — x^XpW^ 
for the resolved case. In these equations we denote the integration regions of the 
direct and resolved convolutions, which are defined by the experimental cuts, by 
and Rf. The parton distribution for a parton of flavour / in the photon and the proton 
are denoted by 0//^(a;^, and (f)f/p{xp,fif) respectively and flnally and 
are the hard scattering functions which, at lowest order, start from the Born cross 
sections. These are the functions that will contain the logarithmic enhancements 
of Qci/Q, and hence depend on the definition of the gap and the gap energy 
flow Qq. We assume that Qq is sufficiently soft that we can ignore the effects of 
emission on the parent jet, known as recoil, but large enough that Qf^ ^ ^qcd- 
The index / denotes the process fj + fp^fi + /2 and the index denotes the 
process 7 + /p — ^ /i + /2- Since the aim of this paper is to calculate ratios of cross 
sections and compare with data, we will take the renormalisation scale to equal the 
factorisation scale and set /if = /ir = Pt, where pt is the transverse momentum of the 
produced jets. 

3.2 Refactorisation 

Following [7, 11] we now refactorise the 2 — > 2 hard scattering function into a hard 
matrix and a soft matrix, 

——{s,t,Arj,as{Hr),Qn,fif) = Hy^{s,t, Arj,as{Hr), fif, fi) 

x4^(Qa,a.(/^.),/^). (3.7) 

We introduce a factorisation scale /i, separate to the parton distribution factorisation 

scale fif, and all dynamics at scales less than n are factored into Slj. Therefore Hjl 
is Qq independent, and all the Qq dependence is included in Sli- This latter function 
describes the soft gluon dynamics. The proof of this statement follows standard 
factorisation arguments [7]. The indices / and L label the basis of colour tensors 
which describe the possible colour exchange in the hard scattering, over which the 
hard and soft matrices are expressed. Soft, wide angle radiation decouples from the 
dynamics of the hard scattering and can be approximated by an effective cross section 
and in this effective theory the partons are treated as recoilless sources of gluonic 
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radiation and replaced by eikonal lines, or path ordered exponentials of the gluon 
field [8]. The soft radiation pattern of this effective eikonal theory then mimics the 
radiation pattern of the partons participating in the hard event, or in other words 
the effective eikonal theory will contain the same logarithms of the soft scale as the 
full theory. The hard scattering function will begin at order for the resolved 
process and order aag for the direct process, and the soft function will begin at 
zeroth order. The lowest order soft function, denoted 5'°^, reduces to a set of colour 
traces. Note that the definition of the gap, and hence the soft function, depends on 
the jet separation Ar/ but we have suppressed this argument of the soft function for 
clarity. 

The construction of the soft function, and in particular its renormalisation prop- 
erties, have been extensively studied elsewhere [7, 8] . A non-local operator is con- 
structed from a product of Wilson lines, which ties four lines (representing the four jet 
process) together with a colour tensor. This operator, which contains ultraviolet di- 
vergences and hence requires renormalisation, is used to construct a so-called eikonal 
cross section, which serves as an effective theory for the soft emission. By summing 
over intermediate states the eikonal cross section is free of potential coUinear sin- 
gularities. It is the ultraviolet renormalisation of the eikonal operator that allows 
colour mixing and the resummation of soft interjet logarithms. 

3.3 Factorisation leads to resummation of soft logarithms 

The partonic cross section, which has been factorised into a hard and a soft function, 
should not depend on the choice of the factorisation scale /x. This leads to the soft 
function obeying 



It is important to point out that we have deliberately ignored the complications of 
terms in this equation arising from radiation into [11], and only include radiation 
emitted by the soft function directly into Vt. The implication of ignoring these non- 
global terms is discussed in section |3.4| , where we also describe how to include their 
effect in a different way. Therefore we have never included the, technically correct, Vt 
argument of the soft function. The matrices Ts{'q,VL) are pro cess- dependent soft 
anomalous dimension matrices that depend on the details of the gap definition and 
the hard scattering. This equation is solved by transforming to a basis in which 
these matrices are diagonal and hence we require a knowledge of the eigenvectors 
and eigenvalues of the soft anomalous dimension matrices. We obtain the entries 
for Ts{'f), VL) from the coefficients of the ultraviolet poles in the matrix of counterterms 
which renormalise the soft function; we can write this quantity as a sum over terms 
from different eikonal lines each with the form of a colour factor multiplied by a 




(3.8) 
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scaleless integral: 

iZs)u = Y.(zl'^hi = J2^tP-^''''- (3-9) 

id id 

The eikonal momentum integrals are process independent, and only depend on i and 
j, the eikonal lines that are connected by the virtual gluon. The colour factor is found 
from consideration of the colour flow for a given process and the basis over which 
the colour flow is to be decomposed. The result is a basis- and process-dependent 
set of colour mixing matrices, which we have listed in appendix P, together with 
our choice of bases in appendix ^ The colour mixing matrices have been obtained 
in [11,23,24] for all relevant subprocesses, and involves using SU(3) colour identities 
like 

t-.tti = ^6,A, - i^f^M (3-10) 

for quark processes and 

da,, = 2 [Tr {m') + Tr {tHh'')] , (3.11) 
fabc = -1i [Tr [tHh') - Tr [m'>)\ , (3.12) 

for gluon processes, to decompose one-loop graphs over a colour basis. We use the 
fact that the colour flow for a real graph is the same as the corresponding virtual 
graph, valid for primary emission. 

Therefore we need to calculate the ultraviolet divergent contribution to the mo- 
mentum function cj^*-'^ from all contributing eikonal graphs. Working in the Feynman 
gauge there are two possible sources of divergence. The first is one loop eikonal graphs 
with a virtual gluon connecting eikonal lines i and j. From the eikonal Feynman 
rules listed in the appendix, these graphs will give a real and an imaginary contri- 
bution to Fg. Note that as we are working in the Feynman gauge the self energy 
diagrams (u;*^**'') give no contribution. The second source of ultraviolet divergences 
are the real emission diagrams, when the emitted gluon is directed out of the gap. 
This can produce an ultraviolet divergence in the eikonal graph as we only measure 
energy flow into the gap and are fully inclusive out of the gap. Hence the virtual 
graphs will only depend on the relative direction of the two eikonal lines and the real 
graphs will give a gap (and hence a jet algorithm) dependence. This sum over real 
and virtual eikonal graphs ensures that the soft function remains free of collinear 
divergences. The imaginary (and geometry independent) part of all our anomalous 
dimension matrices can be extracted from [8,10,11,23], and the calculation for a 
cone-algorithm defined final state has been done in [10]. For the latter case, we have 
re-expressed their results in accordance with our notation in appendix ^ 

By performing the energy integral of the virtual graphs, we can combine the 
result with the corresponding real graph at the integrand level and obtain a partial 
cancellation. We can then write the total momentum part as an integral over the 
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vetoed gap region and arrive at 



+ 5,5,AA-7-^-(l-5,5,), (3.13) 
ZTT ze 

where we integrate over the gap region allowed by the kt algorithm. The notation 5i = 
+ 1(— 1) means that the gluon momentum, denoted k, flows in the same (opposite) 
direction as the momentum of eikonal line i, and Aj = +1 if the eikonal line is a 
quark, or it is a gluon and the soft gluon is above it in the Feynman diagram, or — 1 
if the eikonal line is an antiquark, or it is a gluon and the soft gluon is below it. Pi 
denotes the 4-velocity of eikonal line i, the function 6{k) = 1 when the vector k is 
directed into the gap and in this paper we will use the dimensional regularisation 
convention d = 4 — 2e. The flnite remainder is a result of the energy veto into the 
gap spoiling the real/virtual cancellation. Once we have obtained the momentum 
integrals for the kt deflned flnal state we can construct the anomalous dimension 
matrices using the colour mixing matrices in the appendix. Consideration of the 
eigenvalues and eigenvectors of these matrices, together with the process-dependent 
hard and soft matrices (the full set of hard and soft matrices is shown in appendix 0) 
allows the resummed cross section to be written down, 

^ = E hTsT exp lU^iv, ^) + Mv, ^)) r ^/5o«.(/i)| , (3.14) 



which follows from the diagonalisation of the soft RG equation |3.8| . We denote 
matrices in the diagonal basis by barred matrices, the eigenvalues of the anomalous 
dimension matrices by Aj = a^Aj and we write the lowest-order piece of the QCD beta 
function as f3o = (H^c — 2nj)/(67r). We will observe that, in agreement with [10], 
Re(A) > for all physical channels and hence the resummed cross sections are 
suppressed relative to the fully inclusive cross section. 

3.4 Non-global effects 

As we have discussed in the last section, we have deliberately ignored terms arising 
from secondary radiation into Q, or non-global logarithms (NGLs) [12-15]. Such 
terms arise from radiation at some intermediate scale, M, being emitted outside of Q, 
i.e. into Cl, and then subsequently radiating into Q. In energy flow observables such 
effects give rise to leading logarithms. Inclusion of NGLs in the formalism of the last 
section would result in an explicit M dependence of the soft function and a sensitivity 
to more complicated, 2 ^ n, colour flows for all n > 2. For further details see [11]. 
NG effects have been studied for a two-jet system by Dasgupta and Salam [12,13] and 
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by the current authors with the comphcation of clustering [14] , and in the context of 
energy flow/event shape correlations by Dokshitzer and Marchesini [25] and Berger, 
Kiics and Sterman [26]. 

The effect of NGLs for four-jet kinematics has not been explicitly calculated to 
date. The best that has been managed is a two-jet calculation in the large- A^c limit. 
The NG contributions to the gap cross section factorize into an overall suppression 
factor S'^'^, making it smaller than would be predicted by the resummation of pri- 
mary logarithms alone. In the absence of a complete calculation, we include the 
NGLs approximately, by using our all-order results in the large- iVc limit for S'^'^ in 
a two-jet system [14]. Since four-jet configurations are dominated, in the large- Ai'c 
limit and for large Arj, by colour flows in which two colour dipoles stretch across the 
gap region, we approximate the four-jet NG suppression factor by the square of the 
two-jet one. 

We have reperformed our previous calculation for the kinematic range relevant 
to HERA and find that the variation of S^'-' with Arj is very weak, so we neglect it. 
The variation with is very strong on the other hand. S^'^ is a function of t, 



1 



log 



(3.15) 



where jSo = (IIC^ — 2n/)/(67r), and is well-approximated by a Gaussian in t. Thus 
if is too close to Aqcd, t varies rapidly with it and S'^'^ varies very rapidly. 

It is impossible to quantify the uncertainties in this 
approximation, without a more detailed understanding 
of the underlying physics. To get an idea however, we 
estimate the possible size of higher order corrections, 
by varying the hard scale at which is evaluated. To 
leading logarithmic accuracy, this is equivalent to vary- 
ing the value of as{Q) by an amount proportional to its 



Qn [GcV] 




0.5 




1.0 


47+0-16 

-0.22 


1.5 


fi^+°-^° 

u-oJ-o.ia 


2.0 


'^-0.08 



value. We therefore evaluate t, and hence S (t) , us- 



Table 2: The non-global 
emission suppression factors 
for the 4-jet system, ob- 
tained from an all-orders cal- 
culation for Q = 6 GeV. 



ing our central value of as{Mz) = 0.116, which implies 
as{Q = 6 GeV) = 0.196, and with raised and lowered 
values a^P(6 GeV) = 0.234 and af™'^(6 GeV) = 0.158. 
For Qn = 1.0 GeV, for example, these values result in 

t = Omrto^oll and hence S^^{tf = OA7toH We show the results for all relevant 
values of Qn in table |[ Note that Qn = 0.5 GeV is so low that the range of uncer- 
tainty in if: extends beyond Aqcd and hence the estimate of S^'^ extends to zero. We 
have not shown any results for the 1995 cone-based ZEUS energy flow analysis [4] 
because the low value of Qn = 0.3 GeV means that the central value of the NG 
suppression is already zero, indicating a breakdown of our perturbative approach. 
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The uncertainty on the secondary emission probabihty estimated in this way 
should be added to that on the primary emission probabihty, described in section |[ 
However, we wiU find that the secondary uncertainty generally dominates the two. 
This is therefore clearly an area that needs more work if more precise quantitative 
predictions are to be made. 



4. Soft gluon dynamics for a kt defined final state 

We now evaluate the momentum integral, over the gap region Q. The region 

of integration is determined by the experimental geometry, in which the final state 
is defined by the kt algorithm, and we shall work with the quantity 



where we define k = k/kt- Therefore 

J^^^ = -^A,A,5,5,lnlf + LP.. (4.2) 

We denote the geometry independent imaginary part by LP., and we define the finite 
piece r'^*-'^ by 

^ = (4.3) 
We have extracted the sign function from T^^^\ 

Sij = AiAA6j, (4.4) 

so that 

rfe-) = gfifo-) + J.P.. (4.5) 

In this work we denote the rapidity separation of the jets by At] and the width of an 
azimuthally symmetric rapidity gap by Ay (< Arj). Therefore the available phase 
space for soft gluon emission for a kt defined final state is given by 

= ^lim ^ {^fiAy, Ar^) - n^!'\Ay, A^, R) - d-^\Ay, Ar^, R)) , (4.6) 

where the first term arises from an azimuthally symmetric gap of width Ay, and 
we subtract the region around each jet which is vetoed by the kt algorithm. The 
regions of this equation are shown in figure |l|. In these regions any soft radiation 
is clustered into the jet, and cannot form part of fl. In the first term we take Ay 
approaching Arj, and hence it contains a collinear divergence when the emitted gluon 
is collinear to one of the jets. The two subtracted pieces then remove the regions of 
phase space defined by 

iVk-V^y+i(|)k-(|>^)^<R^ (4.7) 
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f 



Figure 1: The phase space regions for a kt defined final state. The shading denotes the 
regions vetoed by the algorithm, which are subtracted from the piece. Note that we 
have dropped the (ij) superscripts in this figure. 



where the index i labels final state jets and k labels the emitted gluon. The collinear 
divergences in the subtracted pieces exactly match the divergences in the first piece 
and hence the function Q^i^plAt]) is collinear safe. Explicit definitions of the Q 
functions are 

+A./2 ,2, ^_ . p 



(If 



dr] 



(ij) 



-Ay/2 
+Ay/2 

Ar)/2-R 



dr] 



2n{l3,.k){(3,-k) 



(ij) 



(4. 



where we write (/(um = ~ iv ~ A?7/2)2 and obtain ^2 by the symmetry fi, 

where the mapping i is given by {a, b, 1, 2} — > {b, a, 2, 1}. If we define the 
following combinations of momentum integrals, 



where we have combined classes of diagram with the same colour structure, we obtain 
the following closed form for the positive gap contributions, in the limit Ay At], 

/ _ . . OA„. . 1 

TT 

vr 



a 



P 



■ 2Ar/ + log(l - e-^^") + log 



log(l - + log 



7 



TT 



(^-2Ar/ - log(l - e~2Ar,^) _ 



Ar] 


-Ay 




1 


Arj 


-Ay 






Ar] 


-Ay 



2m 



(4.10) 
(4.11) 
(4.12) 



The subtraction pieces are straightforward to express as power series in R and e 
and we shall illustrate the calculation of the momentum integrals with an example. 



4.1 Calculation of ^2^°^^ 

We can write the matrix element in terms of the rapidity of the emitted gluon and 
obtain the following matrix element 



(A • k){Pj ■ k) e-'-i{cosh{Ar]/2 - t]) - cos0) ' 
The integrations for the function Q^"^^ are straightforward, and we obtain 

f ^ ^\smh{Ar]/2- Ay/2)J ^ ' 

The expression for Q^"^^ is 



JAv/2-R J- 



Iav/2-r 7-011^ 27r e ''(cosh(Ar^/2 - r^) - cos 0) ' 

/■+Ay/2 

= / dr]f{r],Ar],R), 

J Ari/2-R 

= r dv'fiv',^V,R), (4.15) 

7Ar?/2-A2//2 

where is defined in the previous section, we have performed the azimuthal in- 
tegration in the second step and changed variable to rj' — Ar]/2 — 77 in the third 
step. The function / can be easily obtained, but it is rather lengthy so we do not 
reproduce it here. We now note that this expression for ^7^*''^ only involves jet 1 and 
hence, by Lorentz invariance, cannot depend on the other jet and so may not be a 
function of the jet separation Ar]. Therefore we write 

^t'^= r dri'f{ri',R). (4.16) 

JAr]/2-Ay/2 

This function /{r]', R) has a divergence as rj' — > 0, so we add and subtract this 
divergence to obtain 



nt'^=f drj'(f{v',R)--)+ f % (4.17) 

JAr]/2-Ay/2 \ V / JAr//2-Aj//2 V 



We can rewrite the lower limit of the first, divergence free, integral as 0, and the 
coUinear divergence is now contained in the second term. Therefore we have used 
Ay as a cut-off for the divergence, and we can write 

^5"^) = + \og2R - log(A77 - Ay). (4.18) 
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We will always denote the divergence free angular integration, which always results 
from such a subtraction, as a barred quantity. We can now rescale the Q^^^^ integral, 
using fj = 7]'/ R, to obtain 

nt''^ = iyv(^R-g{v,R)-^y (4.19) 

This quantity, which is only a function of R, can now be expressed as a power series 
in R and the integrals done on a term-by-term basis. Doing this we obtain the rapidly 
converging series, 

, , , Zlt It It It It It It / . X 

n\ > = -\og(2) + + + + .... 4.20 

^ ^ TT 8 IStt 576 54007r 5292007r 4147200 ^ ^ 

To calculate ^2"^'' we use the parity symmetry and obtain the expression, 

= / dri —— — — ^ r. (4.21) 

Ja,/2-r J-4>r.^ 27re^(cosh(Ar7/2 + r/) + cos0) 

We now perform similar manipulations to the case of However, as Vt^2^^ is a 

function of both final state jets, the resulting expression must be a function of Ar] 
and we also note that ilg"^'* is not divergent. We hence obtain the expression 

= fdr]{R-f{rj,Ar],R)), (4.22) 
Jo 

which we can expand as a power series in the variables R and z = exp(— Ar/), and 
perform the remaining integrations term-by-term. 

The pole arising in the subtraction term ^l^^^^ now cancels against an equivalent 



pole in the function when we expand the latter in Ay around the point Arj, 



lim fij"^^ At] - log( Ar] - Ay) + log(2 sinh Ar/) . (4.23) 

Therefore we find the final, divergence free, form of ^^^^^ as 

Q\^^ = -At] + log(2 sinh Ar]) - log(2i?) - Vl'f^^ - Q^f^ . (4.24) 

We have presented the full set of series expansions in appendix ^ and these, together 
with equations 4.10| - [4-l^ , are sufficient to compute the set of kt defined momentum 



integrals and hence the corresponding anomalous dimension matrix. It is worth 
noting that, although the off-diagonal terms for the kt anomalous dimension matrices 
are no longer pure imaginary as in the cone case, their real parts still vanish for 
large Arj. Indeed for Arj = 2, the real part is more than two orders of magnitude 
smaller than the imaginary part. We have listed the closed- form momentum integrals 
for the cone defined final state using our notation in appendix ^ 
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5. Results 



We now have the tools we need to calculate resummed cross sections at HERA, 
which correctly include primary emission to all orders and secondary emission ap- 
proximately in the large Nc limit. The colour bases used for the contributing partonic 
cross sections are presented in the appendix, along with the decomposed hard and 
soft matrices. We also present the complete colour mixing matrices and the correct 
sign structure for the three classes of diagram. Therefore we can use the eigenvec- 
tors and eigenvalues of the soft anomalous dimension matrices, together with the 
hard and soft matrices, to calculate the primary resummed cross section using equa- 
tion p.l4| , for either a kt or a cone defined final state. The differential cross section, 
in At], can then be computed using the cuts given in section 0, both for the totally 
inclusive cross section (no gap) and for the gap cross section at fixed Qq- The gap 
fraction is then found by dividing the latter quantity by the former. All our results 
are computed using GRV photon parton densities [27] and the MRST proton parton 
densities [28] . We have included an estimate of the theoretical uncertainty in the pri- 
mary resummation by varying the hard scale in the evaluation of a^, while keeping 
the ratio of the hard and soft scales fixed. 



5.1 Totally inclusive ep cross section and the gap cross section 




Arj Arj 

Figure 2: The cross sections for the HI data (left) and the ZEUS data (right), which 
was defined using the kt algorithm with R = 1.0. On both plots the solid line is the 
total inclusive cross section, the dashed line is the gap cross section for Qq = 1 GeV 
with only primary emission included, and the dotted lines indicate the range of theoretical 
uncertainty in the prediction. 



The left hand side of figure |^ shows the totally inclusive dijet cross section for the 
HI analysis and the gap cross section for Qq = 1.0 GeV. We have not shown further 
values of Qq as all the plots show qualitatively the same behaviour. We have cross- 
checked our total inclusive cross section against the Monte Carlo event generator 
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2.0 2.5 3.0 3.5 4.0 2.0 2.5 3.0 3.5 4.0 

A?7 A77 

Figure 3: The gap fractions for the HI analysis with a kt defined final state {R = 1.0), 
at varying Qq. Qq = 0.5, 1.0, 1.5, 2.0 GeV for plots (a), (b), (c) and (d) respectively. 
The HI data is also shown. The solid line includes the effects of primary emission and 
the secondary emission suppression factor. The overall theoretical uncertainty, including 
the primary uncertainty and the secondary uncertainty, is shown by the dotted lines. The 
dashed line indicates the gap fraction obtained by including only primary emission. 



HERWIG [29,30] and we obtained complete agreement for the HI and both the ZEUS 
sets of cuts. In figure |^ the solid curve is the total inclusive cross section, the dashed 
line is the cross section with the primary interjet logarithms resummed and the 
dotted lines show the theoretical uncertainty of the primary resummation, estimated 
by varying as as described above. The inclusion of the primary gap logarithms gives 
a substantial suppression of the cross section; our analysis confirms simple "area of 
phase space" arguments which say that the kt defined final state will have greater 
soft gluon suppression than a cone defined final state due to the increased gap area 
in the [r], 0) plane. This plot for the ZEUS analysis is shown in the right hand side 
of figure ^. 
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2.0 2.5 3.0 3.5 4.0 2.0 2.5 3.0 3.5 4.0 



A?7 A77 

Figure 4: The gap fractions for the ZEUS analysis with a kt defined final state {R = 1.0), 
at varying Qq. Qq = 0.5, 1.0, 1.5, 2.0 GeV for plots (a), (b), (c) and (d) respectively. The 
solid line includes the effects of primary emission and the secondary emission suppression 
factor. The overall theoretical uncertainty, including the primary uncertainty and the 
secondary uncertainty, is shown by the dotted lines. The dashed line indicates the gap 
fraction obtained by including only primary emission. 

5.2 Gap fractions 

The gap fraction is defined as the gap cross section, at fixed Qn, divided by the 
total inclusive cross section. Figure |^ shows the gap fraction for the HI cuts at the 
four experimentally measured values of Qq and figure ^ shows the gap fractions for 
the ZEUS analysis. The solid line is the gap fraction curve obtained by including 
the primary emission and the NG suppression factors of table ^ in the prediction 
for the gap cross section. The dotted lines show the theoretical uncertainty of both 
the primary and secondary emission probabilities, and the dashed line shows the gap 
fraction obtained by including only the primary emission contribution. We find that 
our gap fraction is consistent with the HI values for the measured Qn- The large 
uncertainty in the gap fraction predictions comes from an approximate treatment of 
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the NG suppression and from using perturbation theory at ~1 GeV. Nonetheless, 
this uncertainty is principally in the normalisation of the curves and we expect our 
resummation to accurately describe the shape of the gap fraction curves. 

6. Conclusions 

In this paper we have computed resummed predictions for rapidity gap processes at 
HERA. We include primary logarithms using the soft gluon techniques of CSS, and 
include the effects of NGLs using an overall suppression factor computed from an 
extension of our earlier work. The kt definition of a hadronic final state determines 
the phase space available for soft primary emission and we have computed a set of 
anomalous dimension matrices specific to the geometry of the HI and ZEUS analyses. 
Of course this method can be used for any definition of the gap, provided Q is directed 
away from all hard jets. We then compared our predictions with gaps-between-jets 
data from the HI collaboration and found a consistent agreement. The theoretical 
uncertainty of our predictions is relatively large, and generally dominated by the 
secondary emission uncertainty. However our resummed predictions correctly predict 
the shape of the HI data, and the normalisation agrees within errors. There is a 
suggestion that the dependence is not quite right, with our central Qq — 0.5 GeV 
prediction below data and our central Qa — 2.0 GeV prediction above data, although 
all are within our uncertainty. It is possible that a more complete treatment of 
the perturbative/non-perturbative interface would improve this. We expect that 
calculation of primary emission will be correct if Qq is not too large, so that we can 
neglect jet recoil. However our calculation is of sufficient accuracy in the region of 
phase space probed at HERA. 

Our treatment of the NGLs is very approximate. For a fuller treatment, it is 
necessary to extend the extraction of the suppression factor to beyond the large N^. 
fimit and overcome the inherent disadvantages of the numerical methods used. For 
the current application, consideration of the four jet system is also necessary. We 
reserve the latter extension, in the large Nc limit, for future work. Our calculation has 
not included power corrections [31]. The inclusion of such non-perturbative effects 
is required for a full and correct comparison to the experimental data. Again, we 
reserve this for future work. 

Our calculation involves a numerical integration over all kinematic variables, so 
it would be straightforward to calculate the dependence of the gap fraction on, for ex- 
ample, the fraction of the photon's momentum participating in the hard process, x^. 
This code is available from the authors. 

In conclusion, we have shown that the calculation of primary and secondary 
emission patterns can give a good description of rapidity gap data at HERA. A fuller 
treatment would refine our approximation of NGLs and include power corrections. 
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A. Colour bases 

In this section we present the colour bases used in this work. All the bases in 
this section have appeared elsewhere [10,11,23,24], but we show them here for 
completeness. 

The process qq qq 



The process qq qq 



1 , , 1 



C2 = -^^5al5m + -j5a25hl- (A.2) 



The process qg qg 



Cl = SalSb2, 

C2 = db2c{Tp)la, 

C3 = tfb2c{T^)la. (A.3) 



The processes qq — > gg and gg — > qq 

The process gg — > qq has the basis. 



Cl = SabSl2: 

C2 = dabciTp) 12, 

C3 = tfabc{TF)l2- (A.4) 



To find the basis for qq — > gg, we interchange a 2 and b 1. 
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The process gg — > gg 

The complete basis is 

{ci, C2, C3, Pi, Pgg, Ps^, Pio0 io> -P27} , (A-5) 



where 



Cl = - [/af)c'^12c — '^afec/l2c] 

C2 = - [fabcduc + '^abc/l2c] 

C3 = ^ [/alc'^b2c + '^alc/62c] 

-Pi = o'5al562, 
O 

3 

Pss = -^dalcdh2c) 
Psa ~ '^falcfb2c, 



lOffilO 



- i^ab^U — ^a2^6l) — -falcfb2c, 



1 13 

-P27 = 2 i^ab5l2 + Sa2Sbl) " -5al5b2 " -dalcdb2c- (A. 6) 

The direct processes 

Since there is only one colour structure, these are basis independent. 
B. The hard and soft matrices 

We now show the complete set of hard and soft matrices used in this work. These 
matrices have appeared in a variety of forms in [10, 11,23,24]. In all these equations 
we have set A^c = 3 and have written the coupling scale as /i. Note that all our hard 
matrices differ from the normalisation used in [10,23] by a factor of 7r/(2s) 4tu/s^, 
while they agree with that used in [11]. 

The process qq — ^ qq 

The hard matrix has, in the basis |A.1| , the form 



^(1) ^ 1 / If xi ^X2 \ 

9 s P \ ^X2 X3 / 



where we define 



X2 = 3 



si P 



s2 + u2 U2^^2 

Xs-^^ + ^^^-^-^. (B.2) 

The unequal flavour process qq' qq' is found by dropping the s-channel terms from 
these equations, and the unequal flavour process qq — >• q'q' is found by dropping the 
t-channel terms. The hard matrix for qq qq is found using the transformation 
i ^ u. The corresponding soft matrix for all these processes is 

^<°'-f5.J .V (B.3) 



iW-i) 

The process qq — > qq 

The hard matrix has, in the basis |A.2| , the form 



where we define 



Xi 



t' + P 

To 1 



X2 = 3tt - 



id 



^2 + ^2 i£2^52 



For the process qq' qq' only keep the t-channel terms. The corresponding soft 
matrix is 

The process qg — > qg 



The hard matrix has, in the basis A^, the form 



^(1) _ 1 4tM / i| 

24 2s 



2 I \xi \xi X2 I , (B.7) 



3X2 X2 Xs 



where we define 



Xi = 2 

X2 = 1 



£2 



su 

1 £2 {^2 

2 six s£ 



su It^ , 
;^^3 = 3-4-^ — . B.8 

2su ^ ' 
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The hard matrix for the process qg — > gq is found by the transformation i 
corresponding soft matrix is 



f N,{Nl - 1) 














\K{Ni - 1) y 

The processes qq — > gg and gg — > qq 

In the basis |A.4| the hard matrix for these processes has the form 



H 



A 2s P \ f^' 

5X2 2^^^ 2/^3 , 



where we define 



Xi 

X2 
X3 



tu 



u. The 



(B.9) 



(B.IO) 



The constant A = 9 for the process qq — *• gg and A = 64 for the process 
The matrix for the process gg — > qq is found from the transformation t 
soft matrix is 

^2_i/2iVc^ 



(B.ll) 

99 m- 

-> u. The 
fB.12) 







Nj 



The process — > gg 



The hard matrix, in the basis |A.5| has the block-diagonal form 



^(1) 



OsxS 03x5 



(B.13) 



where the matrix -^5x^5 has the form 



if, 



(1) 1 Q;^(/i)7r Atu 
" 16 2s 1^ 



/ 9X1 

1X2 





9 

ixi 
\xi 

JX2 



3 



9 



X2 -3xi \ 



1X2 



X3 -|X2 





(B.14) 
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and we write 



Xi 

X2 
X3 



1 - 


iu 


si 




¥ ~ 


- — + 


sm' 


si 


iu 


u" 

St 


S2 


V? 




tu 


27 




' su 


liu 



, 1 St 

— -9 Kr + TT^ + TT^ 



2 V st £m 2 su 



(B.15) 



For this process the soft matrix is 



/soooooo o\ 

050000 
5 
1 
8 
8 
000000 20 
000000 27 










(B.16) 



The direct processes 

For both these processes the zeroth order soft factor is unity and the hard functions 



are 



«s«em7r 4:tU j U t 

2s P \i u 



^^ 8 ^O^saemT^ 4:tu f-u, s 
m \-iq[q) ^ gq[q)) = gCg^r ^ ( — + 



(B.17) 



where Cq is the electric charge of quark flavour q, in units of the electron charge. 
Note that if the sum for — > qq is taken to be over four flavours, then this gives a 
factor of 10/9. 



C. Colour decomposition matrices 

We now give the full set of colour decomposition matrices, and the sign function S, 



defined by equation ^^41 , for a, /3 and 7, defined by 



7 = 5,ir(''i)+5.2r('^2). 



(C.l) 
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The process qq qq 



-^(q; + 7) 



The signs are 



a + 7 CpOi - j^{a + P + 2^) 



Sa = +1, 

= +1, 
= -1. 



The process qq qq 



yqq-*qq 



-^(q; + 7) 



The signs are 



a + 7 Cp-f - ^{2a + P + -f) 



= +1, 
= +1. 



The process qg — > qg 



V 





-(a + 7) 



X 



iVg-4 



(a + 7) 



f(« + 7) 
X / 



The signs are 



Sa = +1, 
Sf3 = +1, 

Sj — —1, 



and we define 



X 



(a-7)-^r(»i) + ^r(''2). 



The processes qq — > gig and gg — > 

For qq — > gigi we have 



/ Cj.r(«^) + c^r^i^) + 7) 



X' 



iVg-4 

4Arc 



f (/3 + 7) 



(/9 + 7) X' 



The signs are 



and we define 



Sa = +1, 
Sf3 = +1, 

Sj — —1, 



For gg — > qq we have 



(^99- 



>qq 



\ 





(/9 + 7) 



jV|-4 



f (/3 + 7) 
(/? + 7) X" J 



X" 



The signs are 



and we define 



Sa = +1, 
Sf3 = +1, 

Sj — —1, 

1 



^" = >-^>- 2N. 



p(12) , :(ycp(a6) 

2 



The process gg gg 

where the matrix Alsxs is 

A^3x3 = 



•^3x3 Oaxs 
05x3 AlsxS 



f(Q; + /3) 








f(«-7) 








f(/3-7),. 



and the matrix M^xb is 
/ 3/? 



5x5 







|(a + 7) ' |(a + 7) f(a + 2/3-7) ' 








(C.16) 
(C.17) 
(C.18) 

(C.19) 



(C.20) 



(C.21) 
(C.22) 
(C.23) 

(C.24) 



(C.25) 



(C.26) 



3(q + 7) 

=^ (a + 2/3 -7) |(q; + 7) |(a + 7) 

|(a + 7) 








for Nc = 3. The signs are 



|(« + 7) 




f(a-7) + 

|(q; + 7) |(q; + 7) 2Q;-/?-27y 



(C.27) 



Sa = +1, 


(C.28) 


^5/3 = +1, 


(C.29) 


= -1. 


(C.30) 
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The direct processes 



This processes has no matrix structure. 



2Nc 2 




_^p(62) 

2Nc 2 




(C.31) 



D. The r(^J) 



series expansions 



We have not found a closed form for these integrals, but they are straightforward to 
express as power series in R and e~^'^ (by Lorentz invariance, only the contributions 
from dipoles containing jet 2 are Ajy-dependent) , 



(+0.31831 R + 0.06250 7?^ + 0.00884 R^ + 0.00087 i?^ + 0.00003 R^) + 

(-0.08616 i? - 0.03383 - 0.01197 i?^ - 0.00282 i?^ 

-0.00039 R^ + 0.00001 R'^)e-^^''-^'^ + 

(+0.01166 R + 0.00916 R^ + 0.00551 R^ + 0.00305 i?^ 

+0.00122 R^ + 0.00038 R^ + 0.00011 i?^ + 0.00003 R^)e-^^'^'^-^'^ + 

(-0.00158 R - 0.00186 R^ - 0.00162 R^ - 0.00139 i?^ 

-0.00088 R^ - 0.00041 i?^ - 0.00017 i?^ - 0.00007 i?^ - 0.00002 i^^je'^^^''-^) + 

(+0.00021 R + 0.00034 R^ + 0.00039 R^ + 0.00045 i?^ 

+0.00039 R^ + 0.00024 R^ + 0.00013 i?^ + 0.00007 R^ 

+0.00003 i?^ + 0.00001 i?^°)e-^(^''-2) + 

(-0.00003 R - 0.00006 R^ - 0.00008 R^ - 0.00012 i?^ 

-0.00013 i?^ - 0.00010 i?^ - 0.00007 i?^ - 0.00004 R^ 

-0.00003 i?'^ - 0.00001 i?^°)e-^(^^-^) + 

(+0.00001 R^ + 0.00002 R^ + 0.00003 i?^ + 0.00004 R^ 

+0.00003 R^ + 0.00003 R'^ + 0.00002 R^ + 0.00001 R^)e-^^^''-^A, 



- 0.31831 R + 0.06250 i?^ - 0.00884 R^ + 0.00087 i?^ - 0.00003 R^ 




(D.l) 



(D.2) 




(D.3) 




(D.4) 



(+0.00458 R^ + 0.00389 R^ + 0.00229 i?^ + 0.00104 R^ 
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(a2) 



+0.00038 + 0.00012 + 0.00003 R^)e-'^^^'i-^^ + 
(-0.00124 R^ - 0.00158 i?^ - 0.00124 R'^ - 0.00079 R^ 
-0.00041 i?^ - 0.00018 i?^ - 0.00007 i?^ - 0.00002 R^)e-^^'^'^-^'^ + 
(+0.00025 R^ + 0.00043 R^ + 0.00042 i?^ + 0.00034 R^ 

+0.00024 i?^ + 0.00014 i?^ + 0.00007 R^ + 0.00003 R^ + 0.00001 R^^)e-^^^'^-^^ + 
(-0.00005 R^ - 0.00010 R^ - 0.00011 i?^ - 0.00011 R^ 

-0.00010 i?^ - 0.00007 i?^ - 0.00004 i?^ - 0.00002 R^ - 0.00001 i?^°)e-^('^^-^) + 
(+0.00002 R^ + 0.00003 i?^ + 0.00003 R^ + 0.00003 i?^ 
+0.00003 R^ + 0.00002 R^ + 0.00001 R^)e-^^^''-^^^ , 

(+0.06767 R'^ + 0.02872 i?^ - 0.00096 R^ + 0.00002 i?^)e-(^^-^) + 
(-0.01374 - 0.01166 i?^ - 0.00229 i?^ - 0.00038 R^ 
-0.00021 R'^ - 0.00003 i?^)e-^(^''-2) + 
(+0.00248 R'^ + 0.00316 R^ + 0.00124 i?^ + 0.00032 R^ 
+0.00041 i?^ + 0.00027 i?^ + 0.00007 i?^ + 0.00001 R^)e-^^^''-^'^ + 
(-0.00042 R^ - 0.00071 R^ - 0.00042 i?^ - 0.00019 R^ 

-0.00024 i?^ - 0.00019 i?^ - 0.00007 i?^ - 0.00002 R^ - 0.00001 i?^°)e"^(^''~^^ + 
(+0.00007i?2 + 0.00014 R^ + 0.00011 i?^ + 0.00007 i?^ 

+0.00010 R^ + 0.00009 i?^ + 0.00004 R^ + 0.00002 R^ + 0.00001 R^^)e-^^^'T^'^ 
(-0.00001 R^ - 0.00003 R^ - 0.00003 - 0.00002 R^ 
-0.00003 R^ - 0.00004 R'^ - 0.00002 i?^ - 0.00001 R^)e-^^'^'^-^'^) , 
a 



+ 



= -f--iogi?--iog-— ^ 

^ 7rV2^ 2^A?7-Ay 

- 0.31831 R - 0.06250 R^ - 0.00884 i?^ - 0.00087 i?^ 



(D.6) 
0.00003 i?^) , 



where all coefficients larger tfian 10~^ are shown (recall that we are mainly interested 
in the case R — 1, At] > 2). By symmetry, we have Jlg''^ = where the mapping 
i — > z is given by {a, b, 1, 2} — > {b, a, 2, 1}. 



E. The ^^f^^ angular integrals for a cone geometry 

We present these results as they have not appeared previously in this form. They 
have the expression, 



= / dr] 

J-Ay/2 Jo 



n, ■ :h 



27r{(3i-k){(3j-~ky 



(E.l) 
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where the integrand is found from the appropriate 4-momenta, and the phase space 
is taken to be of width Ay and azimuthally symmetric. Note that these expressions 
do not include the sign factors. We obtain 



References 

[1] G. Marchesini and B.R. Webber, Phys. Rev. D 38 (1988) 3419 

[2] Yu.L. Dokshitzcr, V. Khozc and S. Troyan, in Physics in Collision VI, Proceedings 
of the International Conference, Chicago, Illinois, 1986, edited by M. Derrick (World 
Scientific, Singapore, 1987), p. 365. 

[3] J. D. Bjorken, Phys. Rev. D 47 (1993) 101. 

[4] M. Derrick et al. [ZEUS Collaboration], Phys. Lett. B 369 (1996) 55 
[arXiv:hep-ex/9510012]. 

[5] ZEUS Collaboration, The XXXIth International Conference on High Energy 
Physics, Amsterdam, July 2002 Abstract number 852. 

[6] C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 24 (2002) 517 
[arXiv:hep-ex/0203011]. 

[7] N. Kidonakis, G. Oderda and G. Sterman, Nucl. Phys. B 525 (1998) 299 
[arXiv:hep-ph/9801268]. 

[8] N. Kidonakis, G. Oderda and G. Sterman, Nucl. Phys. B 531 (1998) 365 
[arXiv:hep-ph/9803241]. 

[9] G. Oderda and G. Sterman, Phys. Rev. Lett. 81 (1998) 3591 [arXiv:hep-ph/9806530]. 

[10] G. Oderda, Phys. Rev. D 61 (2000) 014004 [arXiv:hep-ph/9903240]. 

[11] C. F. Berger, T. Kucs and G. Sterman, Phys. Rev. D 65 (2002) 094031 
[arXiv:hep-ph/0110004]. 




(E.2) 



29 



[12] M. Dasgupta and G. P. Salam, Phys. Lett. B 512 (2001) 323 [arXiv:hep-ph/0104277]. 

[13] M. Dasgupta and G. P. Salam, JHEP 0203 (2002) 017 [arXiv:hep-ph/0203009]. 

[14] R. B. Appleby and M. H. Seymour, JHEP 0212 (2002) 063 [arXiv:liep-ph/0211426]. 

[15] R. B. Appleby and G. P. Salam, arXiv:hep-ph/0305232, presented at 38th 

Rencontres de Moriond on QCD and High-Energy Hadronic Interactions, Les Arcs, 
Prance, 22-29 March 2003. 

[16] J. C. Collins, D. E. Soper and G. Sterman, Nucl. Phys. B 308 (1988) 833. 

[17] M. G. Sotiropoulos and G. Sterman, Nucl. Phys. B 419 (1994) 59 
[arXiv:hep-ph/9310279]. 

[18] H. Contopanagos, E. Laenen and G. Sterman, Nucl. Phys. B 484 (1997) 303 
[arXiv:hep-ph/9604313]. 

[19] S. Catani, Yu. L. Dokshitzer, M. H. Seymour and B. R. Webber, Nucl. Phys. B 406 
(1993) 187. 

[20] S. D. Ellis and D. E. Soper, Phys. Rev. D 48 (1993) 3160 [arXiv:hep-ph/9305266]. 

[21] J. M. Butterworth, J. P. Couchman, B. E. Cox and B. M. Waugh, Comput. Phys. 
Commun. 153 (2003) 85 [arXiv:hep-ph/0210022]. 

[22] C. F. Weiszacker, Z. Phys. 88 (1934) 612 ; E. J. WiUiams, Phys. Rev. 45 (1934) 729 

[23] N. Kidonakis and J. F. Owens, Phys. Rev. D 63 (2001) 054019 
[arXiv:hep-ph/0007268]. 

[24] C. F. Berger, SUNY PhD thesis, arXiv:hep-ph/0305076. 

[25] Yu. L. Dokshitzer and G. Marchesini, JHEP 0303 (2003) 040 
[arXiv:hep-ph/0303101]. 

[26] C. F. Berger, T. Kucs and G. Sterman, Phys. Rev. D 68 (2003) 014012 
[arXiv:hep-ph/0303051] . 

[27] M. Gluck, E. Reya and A. Vogt, Phys. Rev. D 46 (1992) 1973. 

[28] A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Phys. Lett. B 443 
(1998) 301 [arXiv:hep-ph/9808371]. 

[29] G. Corcella et al, JHEP 0101 (2001) 010 [arXiv:hep-ph/0011363]. 

[30] G. Corcella et al, "HERWIG 6.5 release note", arXiv:hep-ph/0210213. 

[31] G. P. Korchemsky and G. Sterman, Nucl. Phys. B 555 (1999) 335 
[arXiv:hep-ph/9902341]. 



- 30 - 



